3.48 \(\int \frac {(a+a \sin (e+f x))^3 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^5} \, dx\)

Optimal. Leaf size=77 \[ \frac {a^3 c^3 (A+B) \cos ^7(e+f x)}{9 f (c-c \sin (e+f x))^8}+\frac {a^3 c^2 (A-8 B) \cos ^7(e+f x)}{63 f (c-c \sin (e+f x))^7} \]

[Out]

1/9*a^3*(A+B)*c^3*cos(f*x+e)^7/f/(c-c*sin(f*x+e))^8+1/63*a^3*(A-8*B)*c^2*cos(f*x+e)^7/f/(c-c*sin(f*x+e))^7

________________________________________________________________________________________

Rubi [A]  time = 0.24, antiderivative size = 77, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 3, integrand size = 36, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.083, Rules used = {2967, 2859, 2671} \[ \frac {a^3 c^2 (A-8 B) \cos ^7(e+f x)}{63 f (c-c \sin (e+f x))^7}+\frac {a^3 c^3 (A+B) \cos ^7(e+f x)}{9 f (c-c \sin (e+f x))^8} \]

Antiderivative was successfully verified.

[In]

Int[((a + a*Sin[e + f*x])^3*(A + B*Sin[e + f*x]))/(c - c*Sin[e + f*x])^5,x]

[Out]

(a^3*(A + B)*c^3*Cos[e + f*x]^7)/(9*f*(c - c*Sin[e + f*x])^8) + (a^3*(A - 8*B)*c^2*Cos[e + f*x]^7)/(63*f*(c -
c*Sin[e + f*x])^7)

Rule 2671

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_), x_Symbol] :> Simp[(b*(g*
Cos[e + f*x])^(p + 1)*(a + b*Sin[e + f*x])^m)/(a*f*g*m), x] /; FreeQ[{a, b, e, f, g, m, p}, x] && EqQ[a^2 - b^
2, 0] && EqQ[Simplify[m + p + 1], 0] &&  !ILtQ[p, 0]

Rule 2859

Int[(cos[(e_.) + (f_.)*(x_)]*(g_.))^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((c_.) + (d_.)*sin[(e_.)
+ (f_.)*(x_)]), x_Symbol] :> Simp[((b*c - a*d)*(g*Cos[e + f*x])^(p + 1)*(a + b*Sin[e + f*x])^m)/(a*f*g*(2*m +
p + 1)), x] + Dist[(a*d*m + b*c*(m + p + 1))/(a*b*(2*m + p + 1)), Int[(g*Cos[e + f*x])^p*(a + b*Sin[e + f*x])^
(m + 1), x], x] /; FreeQ[{a, b, c, d, e, f, g, m, p}, x] && EqQ[a^2 - b^2, 0] && (LtQ[m, -1] || ILtQ[Simplify[
m + p], 0]) && NeQ[2*m + p + 1, 0]

Rule 2967

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])*((c_) + (d_.)*sin[(e_
.) + (f_.)*(x_)])^(n_.), x_Symbol] :> Dist[a^m*c^m, Int[Cos[e + f*x]^(2*m)*(c + d*Sin[e + f*x])^(n - m)*(A + B
*Sin[e + f*x]), x], x] /; FreeQ[{a, b, c, d, e, f, A, B, n}, x] && EqQ[b*c + a*d, 0] && EqQ[a^2 - b^2, 0] && I
ntegerQ[m] &&  !(IntegerQ[n] && ((LtQ[m, 0] && GtQ[n, 0]) || LtQ[0, n, m] || LtQ[m, n, 0]))

Rubi steps

\begin {align*} \int \frac {(a+a \sin (e+f x))^3 (A+B \sin (e+f x))}{(c-c \sin (e+f x))^5} \, dx &=\left (a^3 c^3\right ) \int \frac {\cos ^6(e+f x) (A+B \sin (e+f x))}{(c-c \sin (e+f x))^8} \, dx\\ &=\frac {a^3 (A+B) c^3 \cos ^7(e+f x)}{9 f (c-c \sin (e+f x))^8}+\frac {1}{9} \left (a^3 (A-8 B) c^2\right ) \int \frac {\cos ^6(e+f x)}{(c-c \sin (e+f x))^7} \, dx\\ &=\frac {a^3 (A+B) c^3 \cos ^7(e+f x)}{9 f (c-c \sin (e+f x))^8}+\frac {a^3 (A-8 B) c^2 \cos ^7(e+f x)}{63 f (c-c \sin (e+f x))^7}\\ \end {align*}

________________________________________________________________________________________

Mathematica [B]  time = 2.52, size = 283, normalized size = 3.68 \[ -\frac {a^3 (\sin (e+f x)+1)^3 \left (\cos \left (\frac {1}{2} (e+f x)\right )-\sin \left (\frac {1}{2} (e+f x)\right )\right ) \left (315 (A-B) \cos \left (\frac {1}{2} (e+f x)\right )-189 (A-B) \cos \left (\frac {3}{2} (e+f x)\right )+189 A \sin \left (\frac {1}{2} (e+f x)\right )+105 A \sin \left (\frac {3}{2} (e+f x)\right )-27 A \sin \left (\frac {5}{2} (e+f x)\right )-A \sin \left (\frac {9}{2} (e+f x)\right )-63 A \cos \left (\frac {5}{2} (e+f x)\right )+9 A \cos \left (\frac {7}{2} (e+f x)\right )+693 B \sin \left (\frac {1}{2} (e+f x)\right )+483 B \sin \left (\frac {3}{2} (e+f x)\right )-225 B \sin \left (\frac {5}{2} (e+f x)\right )-63 B \sin \left (\frac {7}{2} (e+f x)\right )+8 B \sin \left (\frac {9}{2} (e+f x)\right )+63 B \cos \left (\frac {5}{2} (e+f x)\right )-9 B \cos \left (\frac {7}{2} (e+f x)\right )\right )}{504 c^5 f (\sin (e+f x)-1)^5 \left (\sin \left (\frac {1}{2} (e+f x)\right )+\cos \left (\frac {1}{2} (e+f x)\right )\right )^6} \]

Antiderivative was successfully verified.

[In]

Integrate[((a + a*Sin[e + f*x])^3*(A + B*Sin[e + f*x]))/(c - c*Sin[e + f*x])^5,x]

[Out]

-1/504*(a^3*(Cos[(e + f*x)/2] - Sin[(e + f*x)/2])*(1 + Sin[e + f*x])^3*(315*(A - B)*Cos[(e + f*x)/2] - 189*(A
- B)*Cos[(3*(e + f*x))/2] - 63*A*Cos[(5*(e + f*x))/2] + 63*B*Cos[(5*(e + f*x))/2] + 9*A*Cos[(7*(e + f*x))/2] -
 9*B*Cos[(7*(e + f*x))/2] + 189*A*Sin[(e + f*x)/2] + 693*B*Sin[(e + f*x)/2] + 105*A*Sin[(3*(e + f*x))/2] + 483
*B*Sin[(3*(e + f*x))/2] - 27*A*Sin[(5*(e + f*x))/2] - 225*B*Sin[(5*(e + f*x))/2] - 63*B*Sin[(7*(e + f*x))/2] -
 A*Sin[(9*(e + f*x))/2] + 8*B*Sin[(9*(e + f*x))/2]))/(c^5*f*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^6*(-1 + Sin[
e + f*x])^5)

________________________________________________________________________________________

fricas [B]  time = 0.43, size = 331, normalized size = 4.30 \[ -\frac {{\left (A - 8 \, B\right )} a^{3} \cos \left (f x + e\right )^{5} - {\left (4 \, A + 31 \, B\right )} a^{3} \cos \left (f x + e\right )^{4} + {\left (19 \, A + 37 \, B\right )} a^{3} \cos \left (f x + e\right )^{3} + 4 \, {\left (13 \, A + 22 \, B\right )} a^{3} \cos \left (f x + e\right )^{2} - 28 \, {\left (A + B\right )} a^{3} \cos \left (f x + e\right ) - 56 \, {\left (A + B\right )} a^{3} + {\left ({\left (A - 8 \, B\right )} a^{3} \cos \left (f x + e\right )^{4} + {\left (5 \, A + 23 \, B\right )} a^{3} \cos \left (f x + e\right )^{3} + 12 \, {\left (2 \, A + 5 \, B\right )} a^{3} \cos \left (f x + e\right )^{2} - 28 \, {\left (A + B\right )} a^{3} \cos \left (f x + e\right ) - 56 \, {\left (A + B\right )} a^{3}\right )} \sin \left (f x + e\right )}{63 \, {\left (c^{5} f \cos \left (f x + e\right )^{5} + 5 \, c^{5} f \cos \left (f x + e\right )^{4} - 8 \, c^{5} f \cos \left (f x + e\right )^{3} - 20 \, c^{5} f \cos \left (f x + e\right )^{2} + 8 \, c^{5} f \cos \left (f x + e\right ) + 16 \, c^{5} f - {\left (c^{5} f \cos \left (f x + e\right )^{4} - 4 \, c^{5} f \cos \left (f x + e\right )^{3} - 12 \, c^{5} f \cos \left (f x + e\right )^{2} + 8 \, c^{5} f \cos \left (f x + e\right ) + 16 \, c^{5} f\right )} \sin \left (f x + e\right )\right )}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+a*sin(f*x+e))^3*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^5,x, algorithm="fricas")

[Out]

-1/63*((A - 8*B)*a^3*cos(f*x + e)^5 - (4*A + 31*B)*a^3*cos(f*x + e)^4 + (19*A + 37*B)*a^3*cos(f*x + e)^3 + 4*(
13*A + 22*B)*a^3*cos(f*x + e)^2 - 28*(A + B)*a^3*cos(f*x + e) - 56*(A + B)*a^3 + ((A - 8*B)*a^3*cos(f*x + e)^4
 + (5*A + 23*B)*a^3*cos(f*x + e)^3 + 12*(2*A + 5*B)*a^3*cos(f*x + e)^2 - 28*(A + B)*a^3*cos(f*x + e) - 56*(A +
 B)*a^3)*sin(f*x + e))/(c^5*f*cos(f*x + e)^5 + 5*c^5*f*cos(f*x + e)^4 - 8*c^5*f*cos(f*x + e)^3 - 20*c^5*f*cos(
f*x + e)^2 + 8*c^5*f*cos(f*x + e) + 16*c^5*f - (c^5*f*cos(f*x + e)^4 - 4*c^5*f*cos(f*x + e)^3 - 12*c^5*f*cos(f
*x + e)^2 + 8*c^5*f*cos(f*x + e) + 16*c^5*f)*sin(f*x + e))

________________________________________________________________________________________

giac [B]  time = 0.27, size = 301, normalized size = 3.91 \[ -\frac {2 \, {\left (63 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{8} - 63 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{7} + 63 \, B a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{7} + 483 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{6} + 105 \, B a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{6} - 315 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{5} + 315 \, B a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{5} + 693 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{4} + 189 \, B a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{4} - 189 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{3} + 189 \, B a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{3} + 225 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} + 27 \, B a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )^{2} - 9 \, A a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 9 \, B a^{3} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 8 \, A a^{3} - B a^{3}\right )}}{63 \, c^{5} f {\left (\tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right ) - 1\right )}^{9}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+a*sin(f*x+e))^3*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^5,x, algorithm="giac")

[Out]

-2/63*(63*A*a^3*tan(1/2*f*x + 1/2*e)^8 - 63*A*a^3*tan(1/2*f*x + 1/2*e)^7 + 63*B*a^3*tan(1/2*f*x + 1/2*e)^7 + 4
83*A*a^3*tan(1/2*f*x + 1/2*e)^6 + 105*B*a^3*tan(1/2*f*x + 1/2*e)^6 - 315*A*a^3*tan(1/2*f*x + 1/2*e)^5 + 315*B*
a^3*tan(1/2*f*x + 1/2*e)^5 + 693*A*a^3*tan(1/2*f*x + 1/2*e)^4 + 189*B*a^3*tan(1/2*f*x + 1/2*e)^4 - 189*A*a^3*t
an(1/2*f*x + 1/2*e)^3 + 189*B*a^3*tan(1/2*f*x + 1/2*e)^3 + 225*A*a^3*tan(1/2*f*x + 1/2*e)^2 + 27*B*a^3*tan(1/2
*f*x + 1/2*e)^2 - 9*A*a^3*tan(1/2*f*x + 1/2*e) + 9*B*a^3*tan(1/2*f*x + 1/2*e) + 8*A*a^3 - B*a^3)/(c^5*f*(tan(1
/2*f*x + 1/2*e) - 1)^9)

________________________________________________________________________________________

maple [B]  time = 0.53, size = 205, normalized size = 2.66 \[ \frac {2 a^{3} \left (-\frac {992 A +800 B}{6 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{6}}-\frac {304 A +144 B}{4 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{4}}-\frac {A}{\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1}-\frac {928 A +864 B}{7 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{7}}-\frac {86 A +26 B}{3 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{3}}-\frac {512 A +512 B}{8 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{8}}-\frac {14 A +2 B}{2 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{2}}-\frac {680 A +440 B}{5 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{5}}-\frac {128 A +128 B}{9 \left (\tan \left (\frac {f x}{2}+\frac {e}{2}\right )-1\right )^{9}}\right )}{f \,c^{5}} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((a+a*sin(f*x+e))^3*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^5,x)

[Out]

2/f*a^3/c^5*(-1/6*(992*A+800*B)/(tan(1/2*f*x+1/2*e)-1)^6-1/4*(304*A+144*B)/(tan(1/2*f*x+1/2*e)-1)^4-A/(tan(1/2
*f*x+1/2*e)-1)-1/7*(928*A+864*B)/(tan(1/2*f*x+1/2*e)-1)^7-1/3*(86*A+26*B)/(tan(1/2*f*x+1/2*e)-1)^3-1/8*(512*A+
512*B)/(tan(1/2*f*x+1/2*e)-1)^8-1/2*(14*A+2*B)/(tan(1/2*f*x+1/2*e)-1)^2-1/5*(680*A+440*B)/(tan(1/2*f*x+1/2*e)-
1)^5-1/9*(128*A+128*B)/(tan(1/2*f*x+1/2*e)-1)^9)

________________________________________________________________________________________

maxima [B]  time = 0.54, size = 2701, normalized size = 35.08 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+a*sin(f*x+e))^3*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))^5,x, algorithm="maxima")

[Out]

-2/315*(A*a^3*(432*sin(f*x + e)/(cos(f*x + e) + 1) - 1728*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 3612*sin(f*x +
 e)^3/(cos(f*x + e) + 1)^3 - 5418*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + 5040*sin(f*x + e)^5/(cos(f*x + e) + 1)
^5 - 3360*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 + 1260*sin(f*x + e)^7/(cos(f*x + e) + 1)^7 - 315*sin(f*x + e)^8/
(cos(f*x + e) + 1)^8 - 83)/(c^5 - 9*c^5*sin(f*x + e)/(cos(f*x + e) + 1) + 36*c^5*sin(f*x + e)^2/(cos(f*x + e)
+ 1)^2 - 84*c^5*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 126*c^5*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 - 126*c^5*si
n(f*x + e)^5/(cos(f*x + e) + 1)^5 + 84*c^5*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 - 36*c^5*sin(f*x + e)^7/(cos(f*
x + e) + 1)^7 + 9*c^5*sin(f*x + e)^8/(cos(f*x + e) + 1)^8 - c^5*sin(f*x + e)^9/(cos(f*x + e) + 1)^9) - 15*A*a^
3*(45*sin(f*x + e)/(cos(f*x + e) + 1) - 117*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 273*sin(f*x + e)^3/(cos(f*x
+ e) + 1)^3 - 315*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + 315*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 - 147*sin(f*x
+ e)^6/(cos(f*x + e) + 1)^6 + 63*sin(f*x + e)^7/(cos(f*x + e) + 1)^7 - 5)/(c^5 - 9*c^5*sin(f*x + e)/(cos(f*x +
 e) + 1) + 36*c^5*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 84*c^5*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 126*c^5*s
in(f*x + e)^4/(cos(f*x + e) + 1)^4 - 126*c^5*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 + 84*c^5*sin(f*x + e)^6/(cos(
f*x + e) + 1)^6 - 36*c^5*sin(f*x + e)^7/(cos(f*x + e) + 1)^7 + 9*c^5*sin(f*x + e)^8/(cos(f*x + e) + 1)^8 - c^5
*sin(f*x + e)^9/(cos(f*x + e) + 1)^9) - 5*B*a^3*(45*sin(f*x + e)/(cos(f*x + e) + 1) - 117*sin(f*x + e)^2/(cos(
f*x + e) + 1)^2 + 273*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 - 315*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + 315*sin(
f*x + e)^5/(cos(f*x + e) + 1)^5 - 147*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 + 63*sin(f*x + e)^7/(cos(f*x + e) +
1)^7 - 5)/(c^5 - 9*c^5*sin(f*x + e)/(cos(f*x + e) + 1) + 36*c^5*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 84*c^5*s
in(f*x + e)^3/(cos(f*x + e) + 1)^3 + 126*c^5*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 - 126*c^5*sin(f*x + e)^5/(cos
(f*x + e) + 1)^5 + 84*c^5*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 - 36*c^5*sin(f*x + e)^7/(cos(f*x + e) + 1)^7 + 9
*c^5*sin(f*x + e)^8/(cos(f*x + e) + 1)^8 - c^5*sin(f*x + e)^9/(cos(f*x + e) + 1)^9) - 10*A*a^3*(9*sin(f*x + e)
/(cos(f*x + e) + 1) - 36*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 84*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 - 63*sin
(f*x + e)^4/(cos(f*x + e) + 1)^4 + 63*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 - 1)/(c^5 - 9*c^5*sin(f*x + e)/(cos(
f*x + e) + 1) + 36*c^5*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 84*c^5*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 126*
c^5*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 - 126*c^5*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 + 84*c^5*sin(f*x + e)^6/
(cos(f*x + e) + 1)^6 - 36*c^5*sin(f*x + e)^7/(cos(f*x + e) + 1)^7 + 9*c^5*sin(f*x + e)^8/(cos(f*x + e) + 1)^8
- c^5*sin(f*x + e)^9/(cos(f*x + e) + 1)^9) - 30*B*a^3*(9*sin(f*x + e)/(cos(f*x + e) + 1) - 36*sin(f*x + e)^2/(
cos(f*x + e) + 1)^2 + 84*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 - 63*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + 63*sin
(f*x + e)^5/(cos(f*x + e) + 1)^5 - 1)/(c^5 - 9*c^5*sin(f*x + e)/(cos(f*x + e) + 1) + 36*c^5*sin(f*x + e)^2/(co
s(f*x + e) + 1)^2 - 84*c^5*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 126*c^5*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 -
 126*c^5*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 + 84*c^5*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 - 36*c^5*sin(f*x + e
)^7/(cos(f*x + e) + 1)^7 + 9*c^5*sin(f*x + e)^8/(cos(f*x + e) + 1)^8 - c^5*sin(f*x + e)^9/(cos(f*x + e) + 1)^9
) + 8*B*a^3*(9*sin(f*x + e)/(cos(f*x + e) + 1) - 36*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 + 84*sin(f*x + e)^3/(c
os(f*x + e) + 1)^3 - 126*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 - 1)/(c^5 - 9*c^5*sin(f*x + e)/(cos(f*x + e) + 1)
 + 36*c^5*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 84*c^5*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 126*c^5*sin(f*x +
 e)^4/(cos(f*x + e) + 1)^4 - 126*c^5*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 + 84*c^5*sin(f*x + e)^6/(cos(f*x + e)
 + 1)^6 - 36*c^5*sin(f*x + e)^7/(cos(f*x + e) + 1)^7 + 9*c^5*sin(f*x + e)^8/(cos(f*x + e) + 1)^8 - c^5*sin(f*x
 + e)^9/(cos(f*x + e) + 1)^9) + 42*A*a^3*(9*sin(f*x + e)/(cos(f*x + e) + 1) - 36*sin(f*x + e)^2/(cos(f*x + e)
+ 1)^2 + 54*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 - 81*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + 45*sin(f*x + e)^5/(
cos(f*x + e) + 1)^5 - 30*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 - 1)/(c^5 - 9*c^5*sin(f*x + e)/(cos(f*x + e) + 1)
 + 36*c^5*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 84*c^5*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 126*c^5*sin(f*x +
 e)^4/(cos(f*x + e) + 1)^4 - 126*c^5*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 + 84*c^5*sin(f*x + e)^6/(cos(f*x + e)
 + 1)^6 - 36*c^5*sin(f*x + e)^7/(cos(f*x + e) + 1)^7 + 9*c^5*sin(f*x + e)^8/(cos(f*x + e) + 1)^8 - c^5*sin(f*x
 + e)^9/(cos(f*x + e) + 1)^9) + 42*B*a^3*(9*sin(f*x + e)/(cos(f*x + e) + 1) - 36*sin(f*x + e)^2/(cos(f*x + e)
+ 1)^2 + 54*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 - 81*sin(f*x + e)^4/(cos(f*x + e) + 1)^4 + 45*sin(f*x + e)^5/(
cos(f*x + e) + 1)^5 - 30*sin(f*x + e)^6/(cos(f*x + e) + 1)^6 - 1)/(c^5 - 9*c^5*sin(f*x + e)/(cos(f*x + e) + 1)
 + 36*c^5*sin(f*x + e)^2/(cos(f*x + e) + 1)^2 - 84*c^5*sin(f*x + e)^3/(cos(f*x + e) + 1)^3 + 126*c^5*sin(f*x +
 e)^4/(cos(f*x + e) + 1)^4 - 126*c^5*sin(f*x + e)^5/(cos(f*x + e) + 1)^5 + 84*c^5*sin(f*x + e)^6/(cos(f*x + e)
 + 1)^6 - 36*c^5*sin(f*x + e)^7/(cos(f*x + e) + 1)^7 + 9*c^5*sin(f*x + e)^8/(cos(f*x + e) + 1)^8 - c^5*sin(f*x
 + e)^9/(cos(f*x + e) + 1)^9))/f

________________________________________________________________________________________

mupad [B]  time = 13.39, size = 346, normalized size = 4.49 \[ \frac {2\,\cos \left (\frac {e}{2}+\frac {f\,x}{2}\right )\,\left (\frac {1013\,A\,a^3}{16}+\frac {149\,B\,a^3}{16}-\frac {113\,A\,a^3\,\cos \left (2\,e+2\,f\,x\right )}{4}+\frac {37\,A\,a^3\,\cos \left (3\,e+3\,f\,x\right )}{8}+\frac {7\,A\,a^3\,\cos \left (4\,e+4\,f\,x\right )}{16}-\frac {41\,B\,a^3\,\cos \left (2\,e+2\,f\,x\right )}{4}+\frac {19\,B\,a^3\,\cos \left (3\,e+3\,f\,x\right )}{8}+\frac {7\,B\,a^3\,\cos \left (4\,e+4\,f\,x\right )}{16}+\frac {63\,A\,a^3\,\sin \left (2\,e+2\,f\,x\right )}{8}+\frac {9\,A\,a^3\,\sin \left (3\,e+3\,f\,x\right )}{2}-\frac {9\,A\,a^3\,\sin \left (4\,e+4\,f\,x\right )}{16}-\frac {63\,B\,a^3\,\sin \left (2\,e+2\,f\,x\right )}{8}-\frac {9\,B\,a^3\,\sin \left (3\,e+3\,f\,x\right )}{2}+\frac {9\,B\,a^3\,\sin \left (4\,e+4\,f\,x\right )}{16}-\frac {257\,A\,a^3\,\cos \left (e+f\,x\right )}{8}-\frac {23\,B\,a^3\,\cos \left (e+f\,x\right )}{8}-\frac {63\,A\,a^3\,\sin \left (e+f\,x\right )}{2}+\frac {63\,B\,a^3\,\sin \left (e+f\,x\right )}{2}\right )}{63\,c^5\,f\,\left (\frac {63\,\sqrt {2}\,\cos \left (\frac {e}{2}+\frac {\pi }{4}+\frac {f\,x}{2}\right )}{8}-\frac {21\,\sqrt {2}\,\cos \left (\frac {3\,e}{2}-\frac {\pi }{4}+\frac {3\,f\,x}{2}\right )}{4}-\frac {9\,\sqrt {2}\,\cos \left (\frac {5\,e}{2}+\frac {\pi }{4}+\frac {5\,f\,x}{2}\right )}{4}+\frac {9\,\sqrt {2}\,\cos \left (\frac {7\,e}{2}-\frac {\pi }{4}+\frac {7\,f\,x}{2}\right )}{16}+\frac {\sqrt {2}\,\cos \left (\frac {9\,e}{2}+\frac {\pi }{4}+\frac {9\,f\,x}{2}\right )}{16}\right )} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(((A + B*sin(e + f*x))*(a + a*sin(e + f*x))^3)/(c - c*sin(e + f*x))^5,x)

[Out]

(2*cos(e/2 + (f*x)/2)*((1013*A*a^3)/16 + (149*B*a^3)/16 - (113*A*a^3*cos(2*e + 2*f*x))/4 + (37*A*a^3*cos(3*e +
 3*f*x))/8 + (7*A*a^3*cos(4*e + 4*f*x))/16 - (41*B*a^3*cos(2*e + 2*f*x))/4 + (19*B*a^3*cos(3*e + 3*f*x))/8 + (
7*B*a^3*cos(4*e + 4*f*x))/16 + (63*A*a^3*sin(2*e + 2*f*x))/8 + (9*A*a^3*sin(3*e + 3*f*x))/2 - (9*A*a^3*sin(4*e
 + 4*f*x))/16 - (63*B*a^3*sin(2*e + 2*f*x))/8 - (9*B*a^3*sin(3*e + 3*f*x))/2 + (9*B*a^3*sin(4*e + 4*f*x))/16 -
 (257*A*a^3*cos(e + f*x))/8 - (23*B*a^3*cos(e + f*x))/8 - (63*A*a^3*sin(e + f*x))/2 + (63*B*a^3*sin(e + f*x))/
2))/(63*c^5*f*((63*2^(1/2)*cos(e/2 + pi/4 + (f*x)/2))/8 - (21*2^(1/2)*cos((3*e)/2 - pi/4 + (3*f*x)/2))/4 - (9*
2^(1/2)*cos((5*e)/2 + pi/4 + (5*f*x)/2))/4 + (9*2^(1/2)*cos((7*e)/2 - pi/4 + (7*f*x)/2))/16 + (2^(1/2)*cos((9*
e)/2 + pi/4 + (9*f*x)/2))/16))

________________________________________________________________________________________

sympy [A]  time = 116.40, size = 3262, normalized size = 42.36 \[ \text {result too large to display} \]

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((a+a*sin(f*x+e))**3*(A+B*sin(f*x+e))/(c-c*sin(f*x+e))**5,x)

[Out]

Piecewise((-126*A*a**3*tan(e/2 + f*x/2)**8/(63*c**5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2
268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5
*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/
2 + f*x/2) - 63*c**5*f) + 126*A*a**3*tan(e/2 + f*x/2)**7/(63*c**5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 +
 f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)*
*5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567
*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f) - 966*A*a**3*tan(e/2 + f*x/2)**6/(63*c**5*f*tan(e/2 + f*x/2)**9 - 567*c*
*5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5*f*tan
(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e/2 + f
*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f) + 630*A*a**3*tan(e/2 + f*x/2)**5/(63*c**5*f*tan(e/2 + f*x/
2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)**6 + 7
938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 + f*x/2)**3 - 2268*c**5
*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f) - 1386*A*a**3*tan(e/2 + f*x/2)**4/(63*c**5*f
*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*tan(e/2
+ f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 + f*x/2)
**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f) + 378*A*a**3*tan(e/2 + f*x/2)
**3/(63*c**5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c
**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*ta
n(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f) - 450*A*a**3*ta
n(e/2 + f*x/2)**2/(63*c**5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/
2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 +
5292*c**5*f*tan(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f) +
 18*A*a**3*tan(e/2 + f*x/2)/(63*c**5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(
e/2 + f*x/2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*
x/2)**4 + 5292*c**5*f*tan(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63
*c**5*f) - 16*A*a**3/(63*c**5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f
*x/2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4
 + 5292*c**5*f*tan(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f
) - 126*B*a**3*tan(e/2 + f*x/2)**7/(63*c**5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5
*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e
/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/
2) - 63*c**5*f) - 210*B*a**3*tan(e/2 + f*x/2)**6/(63*c**5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)*
*8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 793
8*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*
tan(e/2 + f*x/2) - 63*c**5*f) - 630*B*a**3*tan(e/2 + f*x/2)**5/(63*c**5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan
(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f
*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2
 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f) - 378*B*a**3*tan(e/2 + f*x/2)**4/(63*c**5*f*tan(e/2 + f*x/2)**9 -
567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5
*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e
/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f) - 378*B*a**3*tan(e/2 + f*x/2)**3/(63*c**5*f*tan(e/2
+ f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*tan(e/2 + f*x/2)*
*6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 + f*x/2)**3 - 226
8*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f) - 54*B*a**3*tan(e/2 + f*x/2)**2/(63*c*
*5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*tan(
e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 + f*
x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f) - 18*B*a**3*tan(e/2 + f*x
/2)/(63*c**5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c
**5*f*tan(e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*ta
n(e/2 + f*x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f) + 2*B*a**3/(63*
c**5*f*tan(e/2 + f*x/2)**9 - 567*c**5*f*tan(e/2 + f*x/2)**8 + 2268*c**5*f*tan(e/2 + f*x/2)**7 - 5292*c**5*f*ta
n(e/2 + f*x/2)**6 + 7938*c**5*f*tan(e/2 + f*x/2)**5 - 7938*c**5*f*tan(e/2 + f*x/2)**4 + 5292*c**5*f*tan(e/2 +
f*x/2)**3 - 2268*c**5*f*tan(e/2 + f*x/2)**2 + 567*c**5*f*tan(e/2 + f*x/2) - 63*c**5*f), Ne(f, 0)), (x*(A + B*s
in(e))*(a*sin(e) + a)**3/(-c*sin(e) + c)**5, True))

________________________________________________________________________________________